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摘要 : 本 文 在 牛顿 动力 学 方程 中 引入 了 逆 耗 散 势 ， 并 研究 孤立 引力 系统 中 物体 的 运动 方程 ， 
发 现在 大 尺度 下 它 可 以 推导 出 与 人 CDM 模型 相似 的 宇宙 膨胀 动力 学 方程 , 以 及 推导 出 盘 状 
星系 的 渐 近 平坦 的 旋转 速度 性 质 , 与 通常 暗物质 模型 不 同 的 是 ,这 种 渐 近 平坦 的 旋转 速度 性 
质 是 时 间 积累 的 结果 ， 而 不 是 力学 直接 作用 的 结果 ， 以 及 讨论 了 Tully-Fisher 关系 ， 发 现 
MOND 模型 中 的 自然 常数 ao 以 及 函数 4 的 形式 在 本 文 的 模型 中 有 着 明确 对 应 的 物理 意义 。 


关键 词 : 暗物质 ， 暗 能 量 ， 牛 顿 动力 学 方程 ， 逆 耗 散 势 ，MOND 


1. 引 言 

暗物质 与 暗 能 量 是 当今 宇宙 学 中 两 大 疑难 [1-4]， 这 两 种 “看 不 见 ” 的 物质 表 
现 出 了 两 种 截然 不 同 的 效应 ， 即 暗物质 产生 引力 ， 而 上 暗 能 量 产 生 反 引力 。 为 了 能 
解释 这 两 种 现象 , 提出 了 各 种 各 样 的 理论 模型 ， 这 些 理论 模型 以 广义 相对 论 为 基 
础 , 主要 包括 为 两 方面 , 一 是 修改 Einstein 场 方 程 的 左边 , 即 修改 时 空 几何 本 身 ， 
这 些 理论 有 (R) 引力 理论 [5-6]、 膜 世界 引力 理论 [7]、MOND[8-9] 等 。 二 是 修改 
Einstein 场 方程 的 右边 ， 即 在 时 空中 加 入 场 或 物质 ， 这 类 理论 有 ACDM [10]、 
Quintessence[11]、phantom[12-13] 等 。 但 是 目前 各 种 理论 仍然 没有 完美 的 解决 所 
有 问题 。 尤其 是 , 关于 暗物质 粒子 的 寻找 实验 和 上 暗 能 量 的 测量 实验 仍然 一 无 所 获 
[14-15]， 这 促使 我 们 想到 ， 也 许 修 改 引 力 方程 是 一 个 可 以 尝试 的 选择 , 一些 学 者 
也 在 这 方面 作 了 很 多 工作 [16-21]。 

在 众多 替代 暗物质 的 模型 中 ， 由 M.Milgrom 所 提出 MOND 模型 属于 较 广泛 
讨论 的 一 个 。 该 模型 提出 了 一 个 经 验 公式 ， 其 中 引入 了 一 个 基本 的 自然 常数 wm， 
能 很 好 的 解释 星系 尺度 的 一 些 观 测 现象 [22-28]， 比 如 解释 Tully-Fisher 关系 等 
[24-25]。 然而 该 模型 仍然 在 解释 其 它 观测 数据 和 理论 自身 上 存在 一 些 问 题 , 例如 
它 不 能 解释 星系 团 中 的 质量 偏差 [29]， 以 及 违反 了 动量 守恒 这 一 基本 的 守恒 定律 
[30]。 总 之 ，MOND 理论 自 提 出 时 并 不 完善 ， 因 而 Bekenstein 等 人 提出 了 相对 论 
性 的 MOND HEW, BH Tensor-Vector-Scalar Theory (TeVeS ) [20,26], T. G. Zlosnik 
等 人 提出 了 Einstein-Aether 理论 [31]， 以 及 M.Milgrom 后 续 又 提出 了 Bimetric 
MOND Theory[32-33]， 这 些 扩展 理论 极 大 的 丰富 了 MOND 的 内 容 ， 但 是 也 存在 
其 他 问题 [34-35]。 

然而 MOND 及 其 扩展 理论 仍然 取得 了 部 分 成 功 ， 这 似乎 暗示 着 它 源 于 一 个 
更 基本 的 理论 体系 。 本文 试图 在 非 相 对 论 牛 顿 动 力学 的 基础 上 , 采用 修正 惯性 的 
理论 方法 [36]， 即 在 牛顿 动力 学 方程 中 引入 耗 散 势 ， 来 讨论 它 在 大 尺度 下 对 物体 
运动 的 影响 ， 并 试图 揭示 MOND 理论 中 自然 常数 a 的 一 些 深层 次 根源 。 


”通讯 作者 : 胡 锦 文 ，Email: xiongbo@whu.edu.cn/2007.hujinwen@163.com. 


2. 带 耗 散 势 能 的 牛顿 动力 学 方程 

在 非 相对 论 情 况 下 ， 和 牛顿 动力 学 方程 为 : 

F=ma (1) 

KP m AY ite, a 为 运动 加 速度 。 

按照 MOND 的 思想 ， 修 改 MOND 理论 通常 有 两 种 途径 : 第 一 种 是 试图 修 
改 方程 的 右边 ， 即 建立 修正 的 惯性 理论 ， 第 二 种 是 修改 方程 的 左边 ， 建立 修正 的 
引力 理论 。 在 相对 论 情况 下 ,惯性 理论 和 引力 理论 是 同等 的 。 本 文 考虑 在 非 相对 
论 情况 下 采用 修正 的 惯性 理论 ， 将 式 (1) 修 改 为 如 下 方程 

F=m(a-Ay) (2) 
其 中 4 为 常数 ，* 为 物体 运动 速度 。 

对 于 一 个 只 有 两 个 物体 组 成 的 孤立 引力 系统 之 ， 其 中 质量 为 六 的 物体 在 质 
量 为 W 的 准 静 止 物体 的 引力 场 中 作 赤 道 平面 上 的 运动 CM >> m )，, 则 根据 式 (2)， 
系统 的 拉 格 朗 日 方程 为 : 

A Gy wey ie 


2 r 
其 中 (x,0) 为 赤道 平面 加 坐标 系 。 
在 式 (2) 中 引入 了 与 速度 相关 的 耗 散 函 数 ， 其 对 应 的 耗 散 势 为 : 


y= -ma +(r0)] (4) 
根据 耗 散 系统 的 哈密 顿 原理 ， 可 得 质量 为 m 的 物体 运动 方程 为 : 
d ,2 Gok 
ae 0)=Ar° 0 


(3) 


2 (5) 

r-rO + sod 

r 

5 JE 4 m 的 物体 沿 着 坐标 系 的 径 问 运动 时 ， 即 96=0， 则 由 式 (5) 的 第 
二 式 可 得 : 


Er 


cir (6) 


r 
RITADE A a 2 FAAET. AREER r= a(t)R CR 
为 宇宙 膨胀 的 随 动 坐标 ，a() 为 宇宙 膨胀 因子 ) ， 并 取 p=poa ”， 以 及 
M =4x13pr 为 半径 为 r 的 球体 内 物质 的 质量 。 则 根据 式 (6) 可 得 : 


r= 


a= — Pa a (7) 
RNA A = ala, MRTE 
a= Po atta (8) 
对 比 ACDM 模型 中 场 方程 的 时 -时 分 量 [10]， 可 得 
~ (9) 


式 (9) 反 映 出 ， 带 北 耗 散 势 的 牛顿 动力 学 方程 在 大 尺度 上 可 以 给 出 与 膨胀 宇 
宙 学 ACDAM 模型 中 相似 的 动力 学 方程 。 
现在 考虑 质量 为 m 的 物体 在 初始 条 件 下 (t=0) 在 质量 为 W 的 准 静 止 物 体 的 


引力 场 中 作 圆 周 运动 ， 即 =0，GM /x? =0，r, 。 在 耗 散 系 


统 中 ，x 、0 会 随 着 


时 间 而 变化 ， 设 + =f()，09=6,g()， 则 式 (5) 中 第 二 式 可 化 为 : 


r- fen 6, = =Ar 
eo GM 1 3 . 
are =A 
7 十 Ca fg°)=Ar 
考虑 在 初始 条 件 下 (= 0) ， 引 力 场 很 弱 时 ， 即 
GM 
二 ~0 
70 
将 式 (12) 代 入 式 (11) 中 ， 则 有 
PEAT 
解 方程 (13) 可 得 
r=nen 


根据 式 (5) 中 的 第 一 式 ， 则 有 : 


27.25 mM 
r 0=r Qe 


根据 式 (14) 和 式 (15)， 可 得 : 
V =V, 
KPV = 97 为 旋转 线 速度 ， 
由 式 (14) 和 式 (160) 可 知 ， 随 着 时 间 的 推移 ,半径 7 
的 物体 旋转 线 速度 却 保持 不 变 。 这 恰好 和 通常 一 些 盘 状 


H 
=l 
FE 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


V, = 0, rp AMIN WEEE 
指数 


旋转 速度 性 质 相 一 致 [38-39]。 
3.MOND 模型 
我 们 知道 ， 在 M.Milgrom 所 提出 的 MOND # 
它 将 牛顿 力学 中 的 实际 加 速度 与 牛顿 加 速度 修改 为 如 下 关系 
ay =au(a/ a) 
其 中 当 ay/a >>1if, u(a/a)~1, a/a, <<1 时 ，A(a 


系 所 呈现 的 渐 近 平坦 的 


论 中 引入 了 一 个 自然 常数 m， 


[8-9]: 
(17) 


/ay)~alayo 


而 从 式 (12) 可 知 ， 当 初始 条 件 的 引力 场 很 弱 时 ， 随 着 时 间 的 推移 ， 质 量 为 六 


的 物体 运动 线 速度 基本 保持 不 变 ， 而 距离 7 却 可 以 仍然 不 断 
对 应 的 强度 为 g。， 则 有 : 


GM 
&o = 2 
70 
将 式 (18) 代 入 式 (16) 中 ， 则 可 得 : 
V*=V; =GMg, 
由 式 (19) 可 以 给 出 : 
Ig(L) = 41g(V) —Ig(Gg, < M/L>) 
其 中 ML 为 质 光 比 。 


式 (20) 正 是 Tully-Fisher 关系 。 对 比 式 (19) 和 MOND 理 


兽 大 。 设 该 弱 引 力 场 


(18) 


(19) 


(20) 


wW, ARI gy =au。 


并 且 由 式 (14) 可 得 物体 运动 的 实际 加 速度 为 : 


21 
r ne” e” e” ( ) 
而 牛顿 加 速度 为 : 
GM GM g a 
ay 二 一 > 二 本 (22) 
r me e 
则 由 式 (21) 和 式 (22) 可 得 
a 1 a 
之 = 二 = 二 (23) 
a e a 


对 比 式 (17) 和 式 (23) 可 见 , 本 文 的 模型 能 够 推导 出 MOND 理论 中 引入 的 经 验 
函数 4 的 形式 。 


4. 结 论 

本 文 尝 试 在 牛顿 动力 学 方程 中 引入 逆 耗 散 势 ， 并 研究 一 个 孤立 的 引力 系统 ， 
发 现 该 修正 后 的 动力 学 方程 可 以 同时 给 出 与 ACDM 模型 相似 的 宇宙 膨胀 动力 学 
方程 和 盘 状 星系 的 近似 平坦 旋转 速度 性 质 ， 而 这 通常 被 认为 是 暗 能 量 和 暗物质 ， 
这 说 明 ， 引 入 了 道 耗 散 势 后 暗 能 量 效应 和 暗物质 效应 是 一 个 物理 本 质 的 两 方面 。 

与 通常 暗物质 模型 不 同 的 是 , 引入 道 耗 散 势 后 盘 状 星系 所 呈现 的 渐 近 平坦 的 
旋转 速度 性 质 是 长 久 时 间 积 累 后 形成 的 结果 ， 而 不 是 直接 的 力学 作用 结果 。 以 及 
在 其 导出 的 宇宙 膨胀 动力 学 方程 中 , 呈现 的 是 一 个 等 效 的 时 变 宇宙 学 常数 ， 而 这 
与 Quintessence 模型 和 phantom 模型 相 一 致 。 

最 后 ， 与 MOND 模型 一 样 ， 该 修正 后 的 牛顿 动力 学 方程 同样 推导 出 了 
Tully-Fisher 关系 , 并 可 以 看 到 MOND 模型 中 引入 的 自然 常数 oa 以 及 函数 4 的 形 
式 在 本 文 的 模型 中 有 一 个 明确 对 应 的 物理 意义 。 但 是 , 本 文 的 讨论 仍然 不 是 一 个 
相对 论 性 的 修正 牛顿 引力 方程 ， 这 还 需 后 续 的 进一步 扩展 研究 。 
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Potential as an Alternative to Dark Matter and Dark Energy 
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Abstract: 

In this paper, the inverse dissipation dispersity is introduced in the Newtonian 
Dynamics equation, and the equation of motion of the object in the isolated 
gravitational system is studied. It is found that at large scale it can derive the 
cosmic expansion dynamics equation similar to the ACDM model and derive the 
asymptotic flat rotational velocity property of spiral galaxies, which it is unlike 
the usual dark matter model that the asymptotic flat rotational velocity property 
is the result of time accumulation rather than the direct effect of mechanics, and 
the Tully-Fisher relationship is discussed and it is found the natural constants 
in the MOND model that a, and the formof the mg function have a clear corresponding 


physical meaning in the model of this paper. 
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